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£j ' Abstract 

CO 

We have developed a series expansion method for calculating the zero-temperature properties of 

s 

\ lattice electron models for variable electron density, i.e. for finite doping away from the half-filled 

q ■ case. This is done by introducing particle fluctuation terms in both the unperturbed Hamiltonian 

and perturbation. The method is demonstrated by application to the 2-chain t — J ladder, where 

[ — ' we provide comparison with previous work as well as obtaining a number of new results. 
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I. INTRODUCTION 



The last 15 years have seen an explosion of interest in the study of lattice models of 
strongly correlated electrons. While, in large part, this has been driven by the field of 
high T c superconductivity, there are other classes of materials whose properties are also 
determined by strong electron correlations, such as Kondo insulators, heavy fermion systems, 
and organic conductors. 

Models for such systems usually start from the well known Hubbard and/or t — J models, 
or generalizations of these. Exact solutions for these models are known only in 1-dimension. 
In higher dimensions, and indeed in 1-dimension for non-solvable variety of analytic 

and numerical techniques have been used: variational methods, exact diagonalizations for 
small systems, quantum Monte Carlo approaches, the density matrix renormalization group 
(DMRG), and perturbative "linked-cluster" series expansions. Yet it is well-known that 
strong-correlated electron models cause difficulties for most numerical techniques. Quan- 
tum Monte Carlo simulations are subject to the infamous "minus sign" problem, while 
exact diagonalization and DMRG calculations suffer from finite-size effects. It is therefore 
important to develop new techniques which may help to shed light on these problem. 

In this paper we focus on series expansions at zero temperatureB. This approach has 
been used to study both the Hubbard mo dei! and the t — J modeffl. In the past this 
method has been restricted to states with one electron per site (half-filling) and to one and 
two-hole excitations away from half-filling. While this region of the phase diagram is itself 
of considerable interest, many of the more interesting aspects of strongly correlated systems 
occur for finite density of holes ("finite doping"). 

We have explored a number of ways in which the series method might be adapted to 
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variable electron density and have been able to obtain consistent results for one model of 
considerable current interest, the t — J laddeiSi'B. The aim of this paper is thus twofold: to 
demonstrate the method, and to present new results for the t—J ladder away from half-filling. 

For completeness we summarize what is known for the t — J ladder, both at half-filling 
and for finite doping. A clear discussion has been given by Troyer, Tsunetsugu and Rice!. 
At half-filling, and in the limit of strong interchain coupling, the ground state consists of 
spin singlet dimers on each rungs. The lowest spin excitation consists of a triplet excitation 
on one rung, propagating via the coupling between rungs. There are also quasiparticle 
excitations (S — |), corresponding in the limit to a hole excitation on a single rung —these 
carry both spin and charge, and separation of spin and charge does not occur in this model. 

The lowest 2-hole excitation consists of a singlet hole pair on one rung, which develop into 
a band of states by propagation along the ladder. The zero-momentum state in this band is 
then the ground state in the 2-hole sector, and then the other states correspond to a gapless 
band of particle-hole charge excitations relative to the ground state. Correspondingly, the 
system is said to be in a C1S0) phase (1 gapless charge mode, gapless spin modes). 

A similar picture appear to hold as one moves away from the strong interchain coupling 
limit, even in the isotropic case which has been the object of most studies. As one moves 
away from half-filling, however, one finds that the lowest triplet excitation corresponds not 
to the simple triplet excitation of a singlet dimer, but to a particle-hole pair excitation on 
different rungs. Thus the triplet gap evolves discontinuously away from half-filling.! 

Poilblanc, Scalapino and Hankei and Miiller and RicJ have presented a plausible phase 
diagram in the n versus J/t plane (for the isotropic case). For J/t % 2 the system is in a 
C1S0 phase (1 gapless charge mode) for low to moderate doping and crosses to a C1S1 phase 
(i.e. the spin-gap vanishes) for higher doping. For J/t > 2 phase separation is predicted to 

3 



occur. For small J/t and low doping Miiller and RiceQ also find ferromagnetic (Nagaoka) 
and C2S2 phases. 

The paper is arranged as follows: in Sect. II we describe our series expansion method; 
Section III and IV present results, for a case with strong interaction coupling and for the 
isotropic case; Section IV gives a summary and discussion. 

II. SERIES EXPANSIONS FOR VARIABLE ELECTRON DENSITY 

The Hamiltonian of the t — J ladder, in usual notation, is 

H — J y^(St,q ■ Si+i )a — -n iia n i+ i i(1 ) + J± y^(Sjj ■ S i)2 — -n^iUi^) 

i,a i 

~t E P(cla,*Ci+l,a,« + H.C.)P -t ± J2 P(4,l,aC h 2,a + H.C.)P (1) 
i,a,a i,a 

where i labels sites along each chain, a (=| or |) and a (=1,2) are spin and leg indices, P is 
a projection operator which excludes doubly occupied sites. J, t are exchange and hopping 
parameters on each chain, while J±, t± are coupling parameters between the two chains, i.e. 
on the rungs of the ladder. 

In our previous work, and also in the present work, we employ a "rung basis", in which 
the 2nd and 4th terms in (|l|) form the unperturbed Hamiltonian, and the remaining terms 
are treated perturbatively. The ground state of the unperturbed Hamiltonian, at half-filling, 
is then a direct product of spin-singlet states on each rung. Spin excitations, at half-filling, 
consist of a spin-triplet on one rung, which propagates coherently along the ladder. One 
and two-hole charge excitations are created by removing one or two electrons from a rung 
state and allowing these to propagate, via the t-hopping term. 

At finite doping the situation is more complex, and it is impossible to identify a suitable 
unperturbed ground state. However it is possible to make progress by relaxing the con- 
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straint on particle number, by adding a particle nonconserving term to the Hamiltonian, 
and by introducing a chemical potential to control the electron density. A variety of particle 
nonconserving terms are possible, but we have found the following form 



fK E ( C i,l,T C i,2,| C i,l,|. C i,2,T + H.c), 
V ^ i a,a 



to give best results. This term creates a spin-singlet on an empty rung or destroys a spin- 
singlet state to create an empty rung. 
Our Hamiltonian, then, is 

H = H + xV (2) 

where 

H = J± Yl( S iA ■ S »,2 - ^n itl n it2 ) ~ t ±J2 P ( C U,a c i,2,a + H.C.)P 

+^ E p (4,iA,2,i - c Iu c Ut + H - c -) p - v EKi + 

V = J E( S V* ' S *+l,a - \n i)a n i+lA ) - t P(cla,*Ci+l,a,a + H.C.)P 
i,a i,a,<x 

-^p(A+A^-^Aa + h ^p (3) 

and x is an expansion parameter. The particle fluctuation term will mix sectors with different 
particle number, but these terms cancel in the physical limit x — 1, and the final results 
should be independent of the field h. In practice h can be adjusted to improve convergence 
for the extrapolation x — > 1. The electron density is determined by the chemical potential 
term, as usual. 

The eigenstates of H are direct products constructed from the nine possible rung states, 
which are given in Table I. For the range of parameters we use, the lowest energy rung state 



is 



|X) = (« - J± - 2/i) | 00) - V2h(\U)- lit)) (4) 



where k is defined in Table 1. This is a spin-singlet state. 

To compute the perturbation series we choose values for the free parameters t, h, /i, and 
derive expansions in powers of x for the quantities of interest. Series have been computed 
to order x 10 for the ground state energy E , electron density n, and the dispersion relations 
for spin-triplet excitations A s (k) and one-hole bonding excitations A C (A;). The calculations 
involve (trivial) 1-dimensional clusters up to 11 rungs, and are limited to this order by 
computer memory constraints. The series are then evaluated at x — 1, corresponding to the 
full Hamiltonian, by integrated differential approximantsi. The series are too numerous to 
be reproduced here, but are available on request. 

III. RESULTS FOR THE DOPED T - J LADDER (STRONG RUNGS) 

As our series expansion is about the rung limit, convergence will naturally be better 
when the rung coupling terms are greater than coupling along the chains. We present in 
this Section results for J±/J — 4. 

We first consider the ground state energy as a function of chemical potential fi. We have 
evaluated series, in powers of x to order x 10 , for various fixed h,fi. The leading terms are 



E /N = --(J ± +2p + K) + x(- - ( J± ± 2 / 2 + Kf J) O(.r') (5) 



with k = y4/i 2 + (Jj_ + 2/i) 2 . Figure [I] shows curves of E versus fi for t = J = 0.25, 
t± = Jj_ = 1. The solid line in this figure is the result for half-filling (n = 1). For this case the 
charge degrees of freedom are frozen out, and the system is equivalent to a Heisenberg spin 
ladder, where the ground state energy per site, for J± = 4 J, is E /NJ± = —0.388037080. 
Including the other terms gives, at half-filling, E Q /NJ = —0.38803708 — 3/16 — //, which 
is the solid line in the figure. For /x > 0.05 our ground state energy is very close to that 



for half-filling. If we take a large value of /x, /x = 5, to make sure the system is half-filled 
at x = 1, the ground state energy is estimated to be E /NJ± = —5.575537(1), which 
agrees with the results for the Heisenberg spin ladder to 6 digits. The results in Figure 
|T] have been obtained by adjusting the field h to obtain best convergence. An illustration 
of the sensitivity of this procedure is shown in Figure [2], where E is plotted versus h, for 
J± — t± — 1, J — t — 0.25, and /x = 0.08. The size of the error bars, which represent the 
spread among different approximants, is least for h ~ 0.25, which is the value used. E itself 
is relatively insensitive to the value of h chosen. 

The electron density is obtained via the standard relation 

n = ~d^-N (6) 
from which we obtain series in x, which are again evaluated at x = 1 via approximants. 
Figure ||] shows curves of electron density versus /x for the same parameters as Fig. [p. Half- 
filling corresponds to large /x, as expected. As /i decreases, the electron density begins to 
drop, and below /x ~ the error bars increase substantially, but the overall trend is very 
clear. 

We have used the data of Figures [1] and | to compute the ground state energy as a 
function of n, and this is shown in Figure In this figure we have subtracted the chemical 
potential term [in from the energy. The ground state energy seems rather insensitive to the 
electron density from these parameters. 

Next we consider the dispersion relation A s (k) for triplet spin excitations. At half- 
filling these excitations correspond physically to a triplet excitation on a rung propagating 
coherently along the ladder. At finite doping the situation is more complicated because 
rung singlets contain an admixture of empty states, and singlet-triplet excitation processes 
include pair creation terms. Nevertheless spin-triplet excitations are well defined and their 



energy dispersion can be calculated via series methodsElO. Figure $ shows the triplet spin 
dispersion curve A s (k) for J± = t± = 1, J = t — 0.25, and /x = 5,0.15,0.08,0.025, which 
correspond to n = 1.000000(1), 0.99(1), 0.90(1), 0.80(2), respectively. For // = 5, the system 
is at half-filling, and we can reproduce the dispersion curve for the Heisenberg spin ladder 
with 4J = J_i_0, up to 6 digits. The figure shows that the triplet dispersion is very sensitive 
to doping, with a significant change already for 1% doping (/i = 0.15). The minimum 
triplet gap remains at k = 7r for low doping (n > 0.8). The size of the gap decreases for 
very low doping, but then begins to increase. By n = 0.8 the dispersion curve has become 
quite flat near the zone boundary, and it appears that the minimum triplet energy shifts to 
intermediate k values for larger doping. 

IV. RESULTS FOR THE DOPED T - J LADDER (ISOTROPIC CASE) 

Most previous studies of the t— J ladder have been for the isotropic exchange case J = J±, 
t = t±. Although our series here are less regular we are able to obtain results which can be 
compared to previous work. We follow, more or less, the order of Section III. 

Figure ^ shows curves of Eq versus /i for t = 0.46, 0.55, 0.75, 1.0. The solid line in this 
figure is the result for half- filling (n = 1) where the ground state energy per site, for J — J±, 
is E /NJ = —0.578043 — 3/8 — /i.0 For \i > 0.1-0.2 our ground state energy is very close to 
that for half-filling, for all t considered. 

Figure [7| shows curves of electron density versus \i for t = 0.55, 0.75, 1, 2 (setting J = 1). 
For n > —0.7 the electron density decreases with increasing t, for constant fi; or equivalently 
for fixed n the value of // increases with t. Around /i ~ —0.8 there is a crossover region 
where the value of n is insensitive to t, and below this value of \x the dependence of n on 
t is the opposite to that discussed above. As t decreases the curves of n versus \i steepen, 
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and at the critical point for phase separation, estimated to be t c = 0.46380, we expect a 
discontinuous drop away from n = 1 to develop at some critical separation value \x c . Our 
series data are not currently expressed as analytical functions of //, and do not allow us to 
explore the phase separation phenomenon any further here. 

We have used the data of Figure [7| to compute the ground state energy as a function 
of n, and these results are shown in Figure ||. For ease of comparison with previous work 
we have subtracted the chemical potential term fin from the energy. Figure § also serves 
to compare our results with previous work. The filled circles at n = 0.5 are our previous 
results! for quarter-filling, for t = 0.55,0.75,1.0. The dashed and solid lines are results 
obtained from a hard-core boson approximation (HCB) and a recurrence-relation method 
(RRM) for t = 2 by Sierra et aZ.EH The crosses are Density Matrix Renormalization Group 
results^, also for t = 2. The latter agree very well with our series estimates, whereas the 
analytic approximations clearly give too high an energy. 

Next we consider the dispersion relation A s (k) for triplet spin excitations. Figure ^ shows 
the triplet spin dispersion curve A s (k) for t = 0.75, and fi = 5, 0.1, -0.1, -0.45, -0.85 which 
correspond to n = 1.000(2), 0.99(3) 0.91(3) 0.61(3), 0.47(4), respectively. The solid curve 
is the dispersion curve for the Heisenberg spin ladder with J = J ± @. This is the dispersion 
curve for half-filling, and agrees very well with the present results for /i = 5. The figure 
again shows that the triplet dispersion is very sensitive to doping, with a significant change 
already for 1% doping (// = 0.1). At small doping the minimum triplet energy remains at 
k = ir, as for half-filling, but for large doping it shifts to intermediate k values. The series 
become more irregular and the error bars are correspondingly larger. For fi = —0.85, the 
shift in the minimum energy has become a dramatic effect. 

Figure 110 shows the triplet spin dispersion curve A s (k) for t = 0.55, 0.75, 1 and 2, and 



fi is chosen to be -0.54, -0.24, 0.14, 1.9 respectively, so that n ~ 0.8. We also show, in the 
Figure, previous results for the minimum energy gap obtained from exact diagonalizations 
by Hayward et As can be seen, the agreement is excellent. 

Finally we turn to charge excitations, in particular the 1-hole bonding excitation (which 
couples rung states 1 and 6, or 1 and 7 in Table 1). These excitations correspond to creation 



or destruction of a hole in the ground state. Figure [TT| show the dispersion A c (k), again for 
t = 0.75 and various /x, together with the corresponding dispersion curve for the Heisenberg 
spin ladder (solid line)@. This curve, which corresponds to half-filling, again agrees very well 
with our present results for \i — 5. A small amount of doping (up to 10%) has little effect 
on the dispersion curve for k < vr/2, but results in a noticeable flattening for larger k. For 
H = —0.45 (corresponding to n ~ 0.6) there is a pronounced minimum at k ~ tt/2 and a 
large increase in excitation energy for larger k. Troyer et al. have previously noted this 
substantial shift in the minimum of the dispersion curve for the bonding and antibonding 
bands, and have identified the minima with the Fermi points of the respective quasiparticle 
bands. 



V. DISCUSSION AND SUMMARY 

The method of linked- cluster perturbation expansions for strongly interacting lattice 
electron models has, until now, been restricted to systems with one electron per site (half- 
filling) or to one or two holes in the half-filling state. In this paper we show, for the first 
time, how this constraint can be relaxed, and we demonstrate its reliability in the case of 
the 2-leg t — J ladder, where comparisons with other methods can be made. The agreement 
is found to be well within the error estimates. 

We present a number of new results for the t — J ladder, for the ground state energy 
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and electron density as functions of chemical potential, and for the ground state energy and 
some excitation energies for different values of electron density. Results are presented for a 
case in which the rung interactions are stronger than the chain interactions (by a factor of 
4) and also for the isotropic case. 

This work opens up new possibilities for exploring other aspects of strongly correlated 
electron models with a finite hole density, and we intend to take up some of these in future 
work. The series approach is not subject to finite-size effects, unlike exact diagonalization, 
the DMRG, or (to a lesser extent) Monte Carlo simulations; and it provides a versatile 
complementary method to these other numerical techniques. It also suffers from some limi- 
tations, of course, for instance, it does not allow one to explore the phase separated regime, 
where the electron density dependence on the chemical potential is discontinuous. 
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FIG. 1: The ground state energy per site Eq/NJ versus fj, for J± = t± = 1, J = t = 0.25. The 
solid line is the result at half-fillino 
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FIG. 2: The ground state energy per site Eq/NJ versus h for J± = t± = 1, J = t = 0.25, and 
H = 0.08. 
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FIG. 3: The electron density n versus fi for J± = t± = 1, J = t = 0.25. 
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FIG. 4: The ground state energy versus electron density n for J± = t± = 1, J = t = 0.25. 
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FIG. 5: The triplet spin dispersion relation for J± = t± = 1, J = t = 0.25 and [i = 
5, 0.15, 0.08, 0.025, corresponding to n = 1, 0.99, 0.9, 0.8. 
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FIG. 6: The ground state energy per site Eq/NJ versus \x for t = 0.4638, 0.55, 0.75 and 1. The 
solid line is the result at half-filling 
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FIG. 7: The electron density n versus \i for t = 0.55, 0.75, 1 and 2. 
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FIG. 8: The ground state energy versus electron density n for t = 0.55 ,0.75, 1 (upper window) and 
2 (lower window). The large solid circles at quarter filling and half- filling are the results of series 
expansions at quarter filling! (for t = 0.55, 0.75 and 1) and series expansions for the Heisenberg 
ladderR The bold long dashed line, solid line and stars are the results of HCB, RRM and DMRGlId 
for t = 2. 
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FIG. 9: The triplet spin dispersion relation for t = 0.75 and different [i. Also presented are results 
for the Heisenberg spin ladder for J = J± (solid line)0. 
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FIG. 10: The triplet spin dispersion relation for t/J = 0.55, 0.75, 1 and 2 and n = 0.8. The full 
solid points are the finite lattice results for the minimum gap [Ref. |i~4| ]. 
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FIG. 11: The dispersion relation for the one-hole bonding excitation at t = 0.75 and different \x. 
Also presented are the results at half-filling (solid line)i. 
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TABLE I: The nine rung states and their energies, where k = ^/4/i 2 + (Jj_ + 2/u) 2 . 



No. 


Eigenstate 


Eigenvalue 


Name 


1 


(k — J\— 2u) I 00) - y/2h(\U)- lit)) 


-(J± + 2fi + k)/2 


lower energy singlet 


2 


III) 

I -4- -4- / 


-2u 

r 


triplet (SL*. = -I) 


3 


i (||i\ + 1 it)) 


-2u 

r 


triplet (SL< = 0) 


4 


ITT) 

III/ 


-2// 

r 


triplet (SL+ = 1) 

1^ V tot / 


5 


+ K + 2a) | 00) - y/2h(\H)- III)) 


(k - J± - 2a) 12 


higher energy singlet 


6 


75(10 1)+ ||0)) 


-t± - n 


electron- hole bonding (S^ ot = —\) 


7 


7l(|0T)+IT0)) 


-t± - n 


electron- hole bonding (S* ot = ^) 


8 


7^(1 o|)- ||0)) 


t±- n 


electron-hole antibonding (S* ot = — §) 


9 


775(1 ot)- IT 0)) 


t±- n 


electron- hole antibonding (S* ot = ^) 



24 



